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{Sj , Abstract: We explore how the colour of any new TeV-scale resonances that decay into 

top quark pairs can be identified by studying the dependence of the observed cross-section 

Q\ | on a central jet veto. To facilitate this study, colour octet resonance production was 

implemented in Pythia8 and colour singlet resonance production is simulated after minor 
modifications. We find that the colour of a 2 TeV resonance can be identified with 10 fb _1 

00 ■ 

• of data at a centre-of-mass energy of 14 TeV for a wide range of couplings, but only if the 

uncertainty in the theoretical prediction is dramatically reduced from its current level. 
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1 Introduction 

The goal of this paper is to illustrate how one can gain insight into the colour of any new 
TeV-scale resonance that might be produced at the LHC. In the case where a resonance 
is observed in the invariant mass spectrum of dijets, the resonance mass will swiftly be 
measured as will its spin. However, it will also be important to identify the colour charge of 
the resonance. Generally speaking, differently coloured resonances will generate a different 
spectrum of accompanying radiation and one should expect to exploit this difference in 
order to establish the colour of the resonance. This matter has been explored in some 
recent papers [1-4] and our purpose is to present a detailed feasibility study that makes 
quantitative statements about the potential of such a measurement at the LHC. 

In this paper, we explore the associated radiation via the introduction of a jet veto. 
In particular, we compute the rate of resonance production subject to the constraint that 
there should be no jets (apart from those arising from the decay of the resonance) lying 
in the central region of rapidity with transverse momentum above some scale, Qq. The 
variation of the cross section with Qq contains important information on the colour of the 
resonance. The method of jet vetoing to probe the colour structure of new physics should 
have other important applications, for example in [5] it was demonstrated that one can 
extract the effective couplings of the Higgs boson to weak vector bosons and to gluons by 
measuring the veto-scale dependence in "Higgs plus dijet" events with a veto on any third 
jet lying between the primary jets. We therefore expect the method to have general utility. 

Jet vetoing has a number of advantages which make it well suited to the task in hand. 
Experimentally, the systematic uncertainties are expected to be small if one focusses on 
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a "gap fraction" (i.e. the ratio of the cross section with a veto to that without). This 
expectation is born out by recent gap fraction measurements in dijet events at the LHC 
[6]. Apart from that, any attempt to identify the colour structure of new physics processes 
will be hindered by background and by pile-up. One can minimize sensitivity to background 
by using observables that are defined over the ensemble of events and which can therefore 
be corrected for background provided it can be subtracted statistically. The gap fraction 
is such an observable since it is a measure of the Qo dependence of the signal cross section. 
Observables that are measured on an event-by-event basis would suffer from inevitable 
contamination from background processes. In the high-luminosity phase of the LHC, there 
will typically be more than 20 pile-up interactions overlaying the signal and any method 
used to extract the colour flow needs to be robust to this extra activity. The jet veto 
method will be robust against pile- up if Qo is chosen to be large and/or the vetoed jets are 
within the acceptance of the inner tracking detectors, as discussed in Section 3.2. 

It is possible to study the colour structure of events using other methods. One such 
possibility is jet scaling [4, 7], which uses the cross-section ratios, a(n + l)/a(n), for specific 
jet multiplicities, n. The information contained in these ratios is very similar to that 
contained in the gap fraction and we do not discuss it any further here. Another possibility 
might be to exploit the structure of jets to probe the colour flow, e.g. the 'jet pull' [2]. 
This approach should yield important information, however it suffers from event-by-event 
contamination from background and the effect of pile-up would also need to be studied. 
We shall focus on the jet veto method in what follows, whilst acknowledging that other 
methods may have a role to play should a new resonance be observed at the LHC. 

Resonant production of a new heavy gauge boson is the primary signal of new physics 
at hadron colliders in several extensions of the Standard Model (SM). One popular scenario 
is within the Randall- Sundrum (RS) model [8-10] where all the SM fields have access to 
the extra dimensional space. In this scenario, the lightest Kaluza-Klein (KK) excitations of 
the gauge bosons primarily decay into top quarks. There are several alternative scenarios 
within this RS framework which have attracted interest as candidates to resolve a number 
of issues with the SM, e.g. the gauge hierarchy problem [8], the fermion mass hierarchy 
[11-13], gauge coupling unification [14-17], providing a dark matter candidate [18, 19]. 
In addition, recent results on the forward-backward asymmetry in top pair events from 
the CDF and D0 experiments indicate a potential deviation from the SM expectation 
[20-23], which could be caused by a colour octet resonance with a mass around 2 TeV 
[24]. For these reasons, we have chosen to investigate how the veto-scale dependence can 
discriminate between a colour octet and a colour singlet resonance, both with a mass of 
2 TeV and spin-1. 1 

The outline of the paper is as follows: In Section 2, we discuss a new implementation 
of heavy gluon resonances in the Pythia8 Monte Carlo event generator. In Section 3, 
we use this implementation to study the feasibility of identifying the colour of such a 
resonance at the LHC. We start from the assumption that a new resonance has been 
observed at the LHC in pp collisions at \fs = 14 TeV. We then show how to quantify the 

x No other colour is possible if the resonance is to couple to quark-anti-quark pairs. 
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probability of sucessfully extracting the colour of the resonance. Finally, we discuss the 
limiting theoretical and experimental uncertainties associated with the measurement. 



2 Implementation of a gluon resonance in Pythia8 

We have implemented the heavy gluon resonance process within the Monte Carlo generator 
Pythia8 [25, 26]. The implementation assumes that the heavy gluon can only be produced 
by quarks, i.e. qq — > G — > qq. At tree-level, the suppression of any coupling to gluons 
is a consequence of the Landau- Yang theorem (for massless gluons and a heavy spin-1 
resonance) [27, 28] and it has also been shown that the suppression persists at one-loop 
in the Bulk RS model [29]. We include interference with SM QCD 2 — > 2 production 
and allow for general couplings to quarks, in particular different couplings to left and 
right-handed states are permitted. Such chiral couplings have recently been discussed in 
connection with the forward-backward asymmetry (Aps) measurements at the Tevatron 
[30, 31]. A forward-backward asymmetry caused by the heavy gluon at leading-order 
requires an axial-vector coupling to both the initial and final state quarks. For this reason, 
such an asymmetry is not expected in the common RS models [32], where the light quark 
couplings are vector-like. However, such studies might be of interest in alternative models. 

Although the implementation is generic, we shall continue to make reference to the 
most common RS models, where the KK gluon coupling to top quarks is significantly 
stronger than to the other flavours and the total width is therefore dominated by the 
width from decays into tops. The fact that the final state flavour is normally different from 
the initial state flavour has the consequence that the light-quark couplings approximately 
determine the production cross section whereas the top coupling fixes the total width. 
The 6-quark coupling is often sufficiently small that it only weakly effects the width and 
cross section, we therefore restrict the parameter space by fixing g\ to either be zero or 
the reference RS model value [9] discussed below. The heavy gluon uses the particle id 
code 5100021 and its corresponding decay channels can be specified through the standard 
Pythia8 particle data scheme. 

Specifically, the total production cross section is given by 

. „. 87rac „ r . , . 

0\S) = 27 S [aSM + OKK + (Tint\ , (2-1) 

ctsm = , (2.2) 
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where the subscripts indicate the contributions from the SM and KK gluon amplitudes and 
their interference. In the above formulae g l J^ and g % P represent the vector and axial-vector 
couplings to the initial and final state quarks, with flavours i and j; the couplings are all 
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relative to the strong coupling in accordance with [33] . The partonic centre-of-mass energy 
is denoted s and the quantities f3j and Aj are defined via 



Pi 




(2.5) 



(2.6) 



The masses and rrij correspond to the KK gluon (which we generally refer to as the 
"heavy gluon") and final state quark masses and the total width is obtained by summing 
the partial widths: 



In addition to specifying the KK gluon mass, the program offers the possibility to 
assign separate values of the KK gluon coupling to light, bottom and top quarks. The most 
common RS models often predict different couplings to different helicity states and therefore 
the couplings can also be assigned separately for left-handed and right-handed quarks. 
Finally, the hard process is integrated within the Pythia8 framework of parton showers, 
underlying event activity, hadronisation and unstable particle decays. The implementation 
has been validated against results in the literature and it has been used to investigate 
characteristic features of a heavy colour octet: some of those results are presented and 
discussed in the appendix. 

3 Extracting the colour of a heavy resonance 

In this section, we investigate the feasibility of distinguishing a 2 TeV colour octet resonance 
from an otherwise identical colour singlet resonance by vetoing on additional jet activity 
outside the di-top system. 

3.1 Simulation and event selection 

Events are generated using the Pythia8 implementation of heavy gluon resonance produc- 
tion, as discussed in Section 2. The couplings used for the gluon resonance are g% = 0.2, 
g\ = and g\ = 3.6. The mass chosen for this study is 2 TeV. The resulting cross section 
x branching-ratio is 1.1 pb and the resonance width is T/M = 0.2. Purely vector couplings 
are assumed and interference with the SM contribution is ignored. Heavy photon produc- 
tion is simulated by changing the colour flow. Specifically, we replace the colour factors 
and a s — >• e^a em , where e q is the electric charge 2 . We also adjusted (i) the coupling of the 
heavy photon to light quarks to reproduce the production rate for the heavy gluon and (ii) 
the coupling of the heavy photon to top quarks in order to match the total width of the 
heavy gluon. The CTEQ5L parton distribution functions are used [34] with the default tune 
to the underlying event (UE) of version 8.130 (Tunel) [25]. Details of the UE tune are 

2 In other words, we assume that the heavy colour singlet couples like a true photon. 




(2.7) 
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Table 1. Parameters of the top-tagging algorithm. Rqa is the distance parameter that is used 
in the Cambridge-Aachen jet finding algorithm. The distance parameter is defined on an event- 
by-event basis since it is dependent on the scalar summed Et of all particles in the event. The 
parameters 5 P and S r are defined in [38]. 

not very important - the sensitivity to non-perturbative physics was tested by repeating 
the analysis, but with hadronisation and multiple parton interactions switched off. We 
observe that the change in the gap fraction is less than 2% and therefore conclude that 
the soft physics modelling is not a crucial component in determining the feasibility of this 
measurement. 

The top quarks will be highly boosted and form a dijet topology and so the primary 
background will come from QCD 2 — > 2 scattering. We define tt and light-jet backgrounds 
separately, the latter being defined with no top-quarks in the final state. In both cases, the 
SM background events were simulated using Pythia8. Top-jet candidates are identified 
using the FastJet library [35-37] and the Johns Hopkins top-tagging algorithm [38] (the 
parameters of the top-tagging algorithm are listed in Table 1). The input to the algorithm 
are stable particles from the MC event record that have \r/\ < 4.9, i.e. within the fiducial 
acceptance of the LHC detectors (neutrinos are removed) . The events are then required to 
contain two tagged top-jets, each of which has pr > 400 GeV, after which the cross section 
for heavy gluon production is 46.8 fb and for heavy photon production it is 31.4 fb. This 
difference demonstrates an interesting and important difference in the tagging efficiency 
for top quarks originating from differently coloured resonances. We have confirmed that 
this difference is not due to non-perturbative physics, in particular we confirm that the 
difference in tagging efficiency arises after the parton shower with hadronization/under lying 
event playing a small role. It follows that, should such a resonance be observed at the LHC, 
the measurement of the colour flow (using a method such as that presented here) will be 
crucial in determining the production cross section. 3 For both resonances, the top tagging 
is sufficient to reduce the background to a manageable level in the region of the resonance. 
Figure 1(a) shows the invariant mass distributions of the tagged top quarks for each type 
of signal and background process. 

In order to compare like-for-like resonances, we make a further adjustment to the heavy 
photon cross section so that it is equal to the heavy gluon cross section after top tagging, 
i.e. do = 46.8 fb. In our final analysis, all results will be presented for a range of values of 
the production cross section (i.e. uq is a baseline value). 

3 The use of jet substructure (e.g. see the review in [39]) should yield additional information regarding 
the colour of the resonance, assuming that the effects of background contamination and pile-up can be 
brought under control. 
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Figure 1. Signal and background distributions after tagging the two leading jets as top candidates, 
(a) The invariant mass of the top candidates, (b) The Ay between the top candidates. 

3.2 Definition of the jet veto region and the gap-fraction 

To identify the additional jet activity, we use the anti-fcx algorithm with R = 0.6. The 
inputs to this second-stage jet finding are the full set of stable interacting particles, i.e. the 
same as to the top-tagging stage. We keep only those jets, j, that are sufficiently far in 
j] — <j) space from the previously tagged top jets, t, i.e. 

AR = (A V l t + A 4*) 1/2 > R ca (3-1) 

where Rqa is defined on an event-by-event basis as discussed in Table 1. A veto can then 
be applied to the remaining jets in order to elucidate the colour structure of the event. 

It might be anticipated that the best separation between the colour singlet and the 
colour octet resonance would occur by vetoing jet activity between the two top-quark jets 
[1]. However, this requires the top-jets to be separated by a large enough rapidity interval. 
Figure 1(b) shows the signal cross sections as a function of the rapidity separation and 
less than 10% of signal events would be retained by a rapidity cut Ay > 2. Fortunately, 
this efficiency loss can be avoided if the veto is applied to jets that lie in a central rapidity 
interval, i.e. \y\ < 1.5. This choice remains sensitive to the different colour structures 
because colour octet resonances radiate preferentially in the forward/backward regions 
whereas colour singlet resonances radiate preferentially in the central region. 

Using a central rapidity interval for vetoing has the added benefit that the jets are 
contained within the acceptance of the tracking systems of the LHC experiments. Even 
though the details of such criteria are beyond the scope of this study, this will be important 
in order to control effects from pile-up, which can be reduced by examining the in-jet track 
activity and determining the amount of this activity that comes from a single vertex. 

The observable of interest is the gap fraction, f gap (Qo), defined as the fraction of the 
events that do not have a jet (in addition to the top-quark decay products) with px > Qo 
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and \y\ < 1.5, i.e. 



WQo)- CT( Q = 3 00 Ge v)- (3 - 2) 

Figure 2(a) shows the gap-fraction as a function of Qo f° r a heavy gluon resonance predicted 
by Pythia8 and also by the matrix element generator MadGraph/MadEvent 5 [40]. The 
Madgraph predictions are determined for up to three partons (in addition to the tt) in 
the final state. Figure 2(b) shows the equivalent curves for a heavy photon resonance. It 
is clear that the Pythia8 and MadGraph predictions differ for each resonance by about 
10% at the lowest values of Qo- Although these theoretical predictions are significantly 
different from each other (we discuss the theory uncertainty further in Section 3.5), Figure 
2(c) shows that the difference between the heavy gluon and heavy photon gap fractions 
remains sizeable in each case and this is all we need for the fidelity of our analysis. We 
therefore use the results obtained using Pythia8 in the following sections to determine 
the sensitivity of the LHC experiments to the colour charge of new resonances. 
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Figure 2. The gap- fraction for heavy gluon (a) and heavy photon (b) resonances predicted by 
Pythia8 and MadGraph. The Madgraph curves correspond to 1,2 or 3 extra partons. (c) The 
ratio of gluon and photon gap-fractions predicted by Pythia8 and MadGraph. 



3.3 Extracting the signal from background 

A pseudo-experiment approach is adopted in order to determine the expected sensitivity to 
the colour of the resonance at a given luminosity, L. For each signal or background process, 
i, the expected number of events is determined by Aj = cr^L, where U{ is the process cross 
section. The actual number of events, rij, that contribute to a given pseudo-experiment is 
determined using a Poisson distribution with mean equal to Aj. The pseudo-experiment 
is constructed by selecting these events at random from reduced MC samples, i.e. after 
top-tagging. 

The tt invariant mass distribution, m t i, is constructed for several values of Qo (i- e - 
after vetoing events that contain additional jets with p T > Qo in the rapidity interval 
\y\ < 1.5). From this we extract the size of the signal as a function of the veto scale for 
20 GeV < Qo < 300 GeV. The number of signal events at each value of Qo is determined 
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by fitting the invariant mass distribution for both signal and background. The signal is 
parameterized by a skewed Breit-Wigner: 

mlT 2 [a + b(m t t - m )] 



(3.3) 



and mo, T, a and b are allowed to vary in the fit. The background is assumed to have 
a shape determined by the tt sample and only its normalization is allowed to vary in the 
fit. Using just the ti background is reasonable since the light-quark background has a 
very similar shape at large invariant masses. Moreover, we have examined in detail what 
happens if we allow the shape of the background fit to vary whilst keeping the shape of 
the true background unchanged: Our final results are very insensitive to even quite large 
changes in this shape and we are in any case confident that an experimental analysis could 
reliably extract the signal cross section from background. After fitting the invariant mass 
distribution, the number of signal events is given by N s = — N^, where Nt is the total 
number of events in the pseudo-experiment that satisfy 1.5 TeV < m t i < 2.5 TeV and Nb 
is the corresponding number of background events in this range, as extracted from the fit. 
The size and location of this mass window would, of course, be optimized in any analysis. 
The gap-fraction as a function of Qq for the signal events can be constructed by 

fgapiQo) = N S (Q = 3W)GeV)' (3 ' 4) 

Figure 3(a) shows the result of a typical pseudo-experiment in the case of a heavy gluon 
signal and assuming an integrated luminosity of 10 fb . The pseudo-data are compared 
to the theoretical predictions (solid lines, Pythia8) for the heavy gluon and heavy photon 
cases. Although there are sizable statistical fluctuations in the pseudo-data 4 , it is clear 
that they are better represented by the heavy gluon prediction. 

Our ability to extract the signal correctly is confirmed in Figure 3(b). Here the gap 
fraction is obtained by computing the mean across all pseudo-experiments. We clearly see 
that the mean gap fractions agree well with the signal-only theoretical predictions for both 
heavy gluon and heavy photon cases. The error bars represent the RMS spread of f gap 
values and indicate that the gap fraction could have significant discriminating power with 
as little as 10 fb -1 of data provided the production rates are not much smaller than we are 
assuming. 

3.4 Extracting the colour 

The colour of the resonance can be obtained using a fit of the form 

fgapiQo) = ai h{Qo) + a 8 /s(Qo) (3.5) 

where f\ (fs) is the pure theory prediction for a colour singlet (octet) resonance and the a, 
are constants that are allowed to vary in the fit to data 5 . Figure 4(a) shows the probability 



4 The errors shown in Figure 3(a) correspond to treating the numerator and denominator in Eq. (3.4) as 
uncorrelated, Poisson, random values. 
J We constrain them to be positive. 
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Figure 3. (a) The data points show the gap fraction in a typical pseudo-experiment for the case of 
a heavy gluon. The curves show the corresponding theoretical predictions for both a heavy photon 
and a heavy gluon; (b) Mean gap fraction obtained after averaging across 1000 pseudo-experiments. 
The errors in (b) represent the RMS spread of values obtained in the pseudo-experiments. 



of obtaining a specific value of as for a heavy gluon signal, P(as\g), assuming an integrated 
luminosity of 10 fb _1 . P(as\g) is strongly peaked at unity, indicating that the gluon is 
identified as such in the majority of pseudo-experiments. Also shown is P(ai|g) and, as 
expected, it is strongly peaked at zero 6 . Figure 4(b) shows the corresponding distributions 
when the true signal is a heavy photon, i.e. P(ai\^/) and P(as\ , y). 

More interesting is the probability that a particular measured value of as is due to a 
gluon resonance and it can be calculated using Bayes' Theorem if we assume that prior to 
this analysis the resonance is equally likely to be a heavy gluon or a heavy photon: 

T3(„\ \ gMg) /q r?\ 

P(gl " 8) = PMg) + P(a 8 |7r (3 - 6) 

Figure 5(a) shows the values of P{g\a%) and P(pf\a&) as a function of as- As expected, 
P(g\ag) increases as as increases, whereas the P(7|as) decreases. Figure 5(b) shows the 
values of P(g|as) for different values of the input luminosity. 

Thus, for each pseudo-experiment a specific value of as can be extracted from the gap 
fraction fit and the probability that the resonance is a gluon (or a photon) can be deter- 
mined using the probability distributions shown in Figure 5: Exactly the same procedure 
could be used with real data. To quantify the feasibility of making such a measurement 
at the LHC, we calculate the fraction of pseudo-experiments that have P(g|as) larger than 
95%, which we denote G95. For example, G95 = 0.77 for the baseline heavy gluon signal 
and 10 fb" 1 of integrated luminosity. It is very likely, therefore, that the LHC experiments 
could identify the colour structure of such a resonance, should it occur in Nature. 



3 We observe that ai + as ~ 1.0 with an RMS variation of less than 2% 
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(a) (b) 

Figure 4. (a) P(<n|g) and P(a$\g) obtained from 10 5 pseudo-experiments for the case of a heavy 
gluon resonance, (b) P(ai|7) and P(as\ r y) obtained from 10 5 pseudo-experiments for the case of a 
heavy photon resonance. Assuming 10 fb _1 of data. 




Figure 5. (a) P(g|a§) and P(7|as) in the case of a heavy gluon resonance and assuming the 
baseline signal size and 10 fb _1 of integrated luminosity, (b) P(g|as) for three different values of 
the luminosity. 

In Figure 6(a) we show how G95 varies as we vary the size of the signal cross section 
relative to the baseline value (erg) and for different integrated luminosities. For signal cross 
sections less than around 1/3 of the baseline value, a luminosity in excess of 50 fb _1 will 
be needed in order to extract the colour the resonance. Figure 6(b) shows the equivalent 
distribution but now for a heavy photon resonance ^95). Figure 7 shows the corresponding 
plots when the probability per pseudo-experiment is increased to 99% (G99 and Tgg) and 
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it is encouraging to note that the plots look rather similar to the previous case. 




Figure 6. The fraction of LHC experiments that would measure P(g\as) > 95% is denoted by 
Gg.5: (a) G95 as a function of the luminosity and signal size assuming a heavy gluon resonance. 
The fraction of LHC experiments that would measure P(-f\ai) > 95% is denoted by r 95 : (b) r 95 
as a function of the luminosity and signal size assuming a heavy photon resonance. The solid line 
in the bottom left corner indicates the region where the significance of any signal is less than 5a 
significance. 




Figure 7. The fraction of LHC experiments that would measure P(g\ag) > 99% is denoted by 
G99: (a) G99 as a function of the luminosity and signal size assuming a heavy gluon resonance. 
The fraction of LHC experiments that would measure P(7|ai) > 99% is denoted by L99: (b) Fgg 
as a function of the luminosity and signal size assuming a heavy photon resonance. The solid line 
in the bottom left corner indicates the region where the significance of any signal is less than 5a 
significance. 
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3.5 Effect of experimental and theoretical uncertainties 

Experimental and theoretical uncertainties will obviously affect our ability to determine 
the colour of the resonance. The experimental uncertainties associated with the extraction 
of the signal will, to a large extent, cancel in the ratio and so will not affect the measure- 
ment in a significant way. These include the top-tagging efficiency, luminosity and energy 
scale/resolution of the top-tagged jets. However, the uncertainties associated with the veto 
jets, i.e. the jet energy scale, jet energy resolution and jet reconstruction efficiency, do af- 
fect our ability to make the measurement. A feel for the likely systematic uncertainty from 
the experiment can be obtained using the ATLAS measurement of dijet production with a 
veto on additional central jet activity [6]. In particular, the total systematic uncertainty on 
the fgapiQo) distributions were at worst 5% for Qq = 20 GeV for dijet systems constructed 
from central high-px jets. 

The current theoretical uncertainties are probably much larger than this. For example, 
the theory predictions in [6] deviate from the ATLAS dijet data by ~ 25% at Qq = 20 GeV 
and in [41] the uncertainty in the theoretical prediction of the gap fraction was found to 
be approaching 50% in some regions of phase space. There is no reason to suppose that jet 
vetoing in boosted ti production is much better understood. Indeed, the difference between 
Madgraph and Pythia8 was shown in Section 3.2 to be about 10%. It is clear, therefore, 
that the theoretical uncertainties are the limiting factor in the feasibility of making this 
measurement. 

We quantify the impact of the theory uncertainty by arbitrarily changing the shape 
of the theoretical predictions used to fit the pseudo-data (i.e. f± and /s). In each pseudo- 
experiment, we choose two uniformly distributed random numbers in the interval [—X, +X], 
which give the fractional shifts that are applied to f\ and f$ at Qq = 20 GeV. The shift 
of each gap fraction at Qq = 300 GeV is by definition zero and for Qq values between 
20 GeV and 300 GeV the shift is obtained by a linear interpolation. Figure 8 shows G95 
for X = 10% and 25% and an uncertainty in the gap fraction at the 25% level clearly 
has a major impact on the measurement. Encouragingly, an overall 10% uncertainty does 
not degrade the measurement much. It is therefore clear that the theoretical uncertainty 
should be reduced to around the 10% level in order to be confident of using jet vetoing as 
a tool to extract the colour of heavy TeV-scale resonances. 
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(a) (b) 

Figure 8. (a) The impact of a 10% uncertainty in the gap fraction shape on the value of G95. (b) 
The impact of a 25% uncertainty in the gap fraction shape on the value of G95. 
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A Appendix 

This appendix highlights some of the features of the heavy gluon implementation in P YTHIA8 
In particular we take a look at the strong enhancement of the ti invariant mass distribution 
at low masses due to parton distribution function (PDF) effects, the potentially large inter- 
ference effects that can occur between the signal and background and the forward-backward 
asymmetry that can be induced by non-chirally symmetry couplings of the heavy gluon to 
quarks. 

For this appendix only, we take as our reference the typical RS couplings: g q L = g g R = 
g b R = —0.2, g\ = g L = 1, g R = 4 [9]. These contain a significant axial- vector coupling to 
tops 7 . In the following, we do also vary the couplings about these reference values. Also 

7 We define g v = (g L + gn)/2 and g a = (ffz - Sh)/2. 
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for this appendix only, when we speak of the tt mass we refer to the invariant mass of the 
top pair at parton level, i.e. before showering and hadronisation. 




Figure 9. (a) Invariant ti mass distribution for a broad gluon resonance, neglecting the SM 
contribution, (b) Comparison of ti invariant mass distributions for several RS-like scenarios. Con- 
tributions from heavy and SM gluon production are included, both from qq annihilation and gluon 
fusion processes. The details are explained in the text. 

The first point of note is that the shape of a wide resonance can be very strongly 
distorted by the PDFs of the incoming particles. The fast increase of the PDFs with 
decreasing x competes with the suppression from the heavy gluon virtuality, enhancing 
the low mass tail of the resonance and suppressing it at high masses. This is illustrated 
in Figure 9(a), which shows the ii mass spectrum for the 'pure' heavy gluon process (i.e. 
without any SM background or interference effects) for a wide resonance with txIq = 3 TeV, 
g\ = 4 and g l a = 0. The shift of the peak can be significant for large resonance mass and 
couplings, but tends not to be larger than about 100 GeV for typical RS scenarios that 
might be seen at the LHC. Similarly, large mass and couplings also imply that a large 
fraction of the resonance cross section resides in the low mass tail and this fraction is 
sensitive to the choice of mass and couplings. It should also be noted that the tail fraction 
can be significant even for a relatively light resonance and that it is quite possible to 
encounter scenarios where only the minority of the cross section resides in the resonance 
peak. Figure 9(a) also demonstrates the small difference introduced by adding some axial 
coupling while maintaining the same overall normalisation (g* = 3.45, g\ = 2.05). It can 
be seen that while the balance of axial and vector couplings does have some effect on the 
low mass tail of the ti mass spectrum, its effect on the main part of the resonance peak is 
very small. 

While most of the SM ti events at the LHC are initiated by gluon fusion, a significant 
fraction comes from quark annihilation, which can interfere with the heavy gluon amplitude. 
The nature of this interference is sensitive to the model parameters and so, in principle, 
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Figure 10. (a) Comparison of invariant tt mass distributions for RS-like and vector coupling 
scenarios. Contributions from both heavy and SM gluon production are included, but only from qq 
annihilation, (b) Forward-backward asymmetry as a function of the tt invariant mass. The details 
are explained in the text. 

could be used to constrain the heavy gluon. Note that the interference changes sign at the 
G mass: RS models typically predict a negative coupling to the incoming light quarks and a 
positive coupling to the outgoing top quarks, which would result in constructive interference 
below the resonance peak and destructive interference above. This would be reversed if 
the couplings to the incoming and outgoing quarks were of the same sign. Figure 9(b) 
compares the tt mass spectrum for the typical RS scenario (g qg < 0) to an equivalent 
scenario with only positive couplings (g qq > 0) and also to the situation where interference 
effects are not taken into account ("No int"). As pointed out, for example in [42], this 
shows that interference effects can certainly be significant, even on top of the full SM tt 
background. We can say a little more about how the strength of the interference depends 
on the balance of axial and vector top couplings: Since only the vector part of the heavy 
gluon coupling interferes with the SM background it is to be expected that reducing the 
vector coupling and increasing axial coupling will reduce the strength of the interference. 
Figure 10(a), which compares the RS scenario with interference that was used to produce 
Figure 9(b) (g* = 2.5, g\ = 1.5) to its purely vector coupling equivalent (g* = 2.92, = 0), 
demonstrates that this difference is small when comparing RS-like to vector scenarios and 
it will be quite negligible once the gluon fusion background is added in (it is not included 
in the figure for clarity). 

To conclude this appendix, we take a look at the forward-backward asymmetry for 
pp collisions at \fs =1.96 TeV. Here we can make direct comparisons with results in the 
literature in order to validate our implementation. Two different sets of parameters were 
considered: 

• mg = 2.75 TeV, g\ = 4.5, g q R = -0.2, g b R = 2, g\ = g\ = and g R = 6.1 based on 
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[30]; 



• tuq = 1.0 TeV, g% = g^=g*=0 and gl = -g\ = -g\ = 1 based on [31]. 

For the heavy gluon contribution to the asymmetry, we find Afb = 0.05 and Afb = 0.14 
respectively (at parton level), which agree with the published results. Figure 10(b) shows 
Afb as a function of the invariant mass of the top pair, using the parameter set from [30]. 
We note that the asymmetry correctly passes through zero at m t i = mg/^/l + 2g%g\. 
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